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We apply the adiabatic self-consistent collective coordinate (ASCC) method to the multi- 
0(4) model and study collective mass (inertia function) of the many-body tunneling motion. 
Comparing results with those of the exact diagonalization, we show that the ASCC method 
succeeds in describing gradual change of excitation spectra from an anharmonic vibration 
about the spherical shape to a doublet pattern associated with a deformed double-well po- 
tential possessing the oblate-prolate symmetry. The collective mass is significantly increased 
by the quadrupole-pairing contribution to time-odd components of the moving mean field. 
In contrast, the cranking (Inglis-Belyaev) mass based on the constrained mean field, which 
ignores the time-odd components, is smaller than the ASCC mass and fails to reproduce the 
exact spectra. 

§1. Introduction 



Q ' Microscopic theories of large-amplitude collective dynamics are a long-standing 

2 ■ and fundamental subject of nuclear structure physics.^ ''^''^ Though many theories 

have been proposed and tested for a variety of phenomena of the large-amplitude col- 
lective motion,^^'-^'-^'-'^'-^'-^'-^"'-^^'-^^'-^^'-^'^'-^^'-^^'-^'^'-^^'-^^'-^"'-^^'-^^'-^^'-^'^'-^^'-^^'-^^-^ 
^ . still many theoretical problems remains not being fully clarified .'^'E^'ES) In par- 

^ I ticular, the microscopic determination of collective mass (inertia of the collective 

motion) is being a difficult and controversial issue. 

The Inglis-Belyaev cranking mass,'!' which is derived by using the adiabatic per- 
turbation theory, has been widely used in the literature. It has been known that 
values of the cranking mass are systematically too low to reproduce experimental 
data of the low-frequency K'^ = 0"*" vibrational modes in deformed nuclei.!^ The 
Inglis-Belyaev cranking mass does not take account of contributions of the residual 
interactions. Especially, effects of the time-odd components of the moving mean- 
field are completely neglected. Importance of such time-odd contributions to the 
collective mass has been discussed in connection with the Thouless-Valatin moment 
of inertia for rotational motions'^'' ''^^ ' and in the context of the adiabatic time- 
dependent Hartree-Fock (ATDHF) theory .ElBED, EDO 

It is also well known that, 

without taking into account the time-odd components, one cannot obtain the cor- 
rect mass for the center of mass motion of a nucleus. In Ref. II ip. the contribution 
from the time-odd components were evaluated on the basis of the ATDHF theory 
and its importance was demonstrated for the isoscalar giant quadrupole modes. In 
this work, however, the pairing correlations were not taken into account. It should 
be noted that, in the time-dependent Hartree-Fock-Bogoliubov (TDHFB) theory 
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for nuclei with superfluidity, the pairing correlations also cause the time-odd com- 
ponents in the moving mean field. To the best of our knowledge, however, there 
are very few paper^^ trying to evaluate the time-odd effects of the pairing correla- 
tions on large-amplitude collective dynamics. This is rather surprising because we 
know that the pairing correlations play crucial roles in low-energy nuclear collective 
dynamics (see e.g., Refs. li2|) . E3|) ). Thus, in this paper, we develop a microscopic 
theory of large-amplitude collective motions in nuclei with superfluidity and dis- 
cuss the time-odd mean-field effects, associated with the pairing correlations, on the 
inertia mass of collective motion. 

We approach to this goal on the basis of the self-consistent collective coordi- 
nate (SCC) methocP^ and its extension to include the pairing correlations.!^ This 
method is based on the TDHFB theory and enables us to extract, in a fully self- 
consistent manner, the optimum collective coordinate and the collective momentum 
from the huge-dimensional TDHFB phase space. This feature of the SCC method is 
in marked contrast to the widely used generator coordinate methoci^ in which the 
collective coordinate (generator coordinate) is chosen in a phenomenological manner. 
The SCC method has been successfully applied to various kinds of low-frequency an- 
harmonic vibrations and high-spin rotational motions.l^J^J^'E3J^'I^J^J^'E3J^'l^^ 
On the other hand, for genuine large-amplitude cases, such as the nuclear fission and 
the shape coexistence phenomena,^^''^^''-^ a practical scheme of solving its basic 
equations had not been available. A possible solution was proposed in Ref. 15^ : 
the adiabatic approximation of the SCC method, called the adiabatic SCC (ASCC) 
method. The ASCC method was first applied to the solvable multi-0(4) model and 
its feasibility was tested.^'^' The multi-0(4) model may be regarded simpli- 
fied version of the pairing plus quadrupole (P-|-Q) interaction model,'^''^''^ and 
it has been used as a testing ground of microscopic theories of nuclear collective 
motions.l^''^''^J^ Recently, Kobayasi et al applied the ASCC method to the 
oblate-prolate shape coexistence phenomena in ^^Se and 72j^],E3,inHJ and successfully 
extracted a collective path connecting the oblate and prolate equilibrium points.!^ 
It turned out that the collective path runs through the triaxially deformed region 
and it almost coincides with the valley line in the deformation-energy surface ob- 
tained by the constrained Hartree-Bogoliubov (HE) calculation. Similar studies have 
been also carried out by Almehed and Walet.'2^''23 In these works, the P-l-Q inter- 
action Hamiltonian was adopted. In order to study the time-odd contributions to 
the collective mass, however, we need to use a more general Hamiltonian including, 
for instance, the quadrupole-pairing interaction, because the time-odd contributions 
from the P-|-Q interactions are known to 

vanish.EJdS 

The importance of the 

quadrupole-pairing interaction for low-frequency collective excitations has been well 
known (see references cited in the review I^H))). 

In this paper, we extend the multi-0(4) model to include the quadrupole- type 
pairing interaction. Varying its strength within a range consistent with experimental 
data for the ratio A2/AQ of the quadrupole and the monopole pairing gaps,'^ we 
evaluate its effect on the collective mass through the time-odd part of the mean 
field. 

This paper is organized as follows: In §2, the basic equations of the ASCC 
method are briefly reviewed. In §3, we apply the ASCC method to a new version 
of the multi-0(4) model with the quadrupole- type pairing interaction. A numerical 
algorithm to solve the ASCC equations for the multi-0(4) model is presented. In §4, 
we present results of numerical analysis. Conclusions are given in §5. A preliminary 
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version of this work was reported previously in Ref . I70|l . 

§2. Basic equations of the ASCC method 

The sec method enables us to extract the "optimal" collective subspace from 
the huge-dimensional TDHFB spaceP^ The large amplitude collective motion 
takes place in this collective subspace. In the case that the collective subspace 
is parametrized by a pair of collective variables, a single collective coordinate q and 
a momentum p, it is called the collective path. In superconducting nuclei, we need 
an additional set of collective variables to take into account the pairing rotational de- 
grees of freedom, which recovers the particle number symmetry broken by the mean 
field approximation.!^ The pairing rotation is described by the particle-number 
variable N and its canonically conjugate gauge angle (p. Thus, a TDHFB state vec- 
tor on the collective path is parametrized by four collective variables {q,p, 93, N). 
In the sec method, the large-amplitude collective motion is assumed to be approx- 
imately decoupled with other non-collective degrees of freedom. This assumption is 
called "the maximal decoupling condition of the collective submanifold." The time 
dependence of the TDHFB state on the collective path is then determined by this 
set of collective variables. 

The time-dependent variational principle for the collective motion is described 

as 

6{4>{q,p,ip,N)\i-^^-H\(^iq,p,^,N)) =0, (2-1) 

where \(p{q,p, ^p, N)) represents a TDHFB state vector on the collective path and the 
variation is taken for all possible deviations around it. The intrinsic state \(l){q,p, N)) 
in the gauge space (associated with the pairing rotation) is defined by 

\cf>{q,p,^,N)) = e-*'^^ mq,P,N)) . (2-2) 

These four collective variables are required to satisfy the following canonical variable 
conditions, 

{(^{q,p,N)\i—\<P{q,p,N))=p, {cP{q,p, N)\ — \<P{q,p, N)) = 0, (2-3) 

{(^{q,p,N)\N\^{q,p,N))=N^No + n, {(^{q,P, N)\ \cl>{q,p, N)) = 0. (2-4) 

Here, A^o denotes the number of particles of the system (we assume a single kind of 
fermions, for simplicity), and the difference n = N — Nq represents the dynamical 
number fluctuation. 

We now introduce the adiabatic approximation to the TDHFB state \(f){q,p, N)): 
We assume that the collective momentum and the number fluctuation are small and 
that the moving state \(j){q,p, N)) is close to \4>{q)) = \4>{q,p = 0,N = Nq)). Using 
the fact that an arbitrary TDHFB state vector can be constructed from a given TD- 
HFB state vector by a unitary transformation (generalized Thouless theorem) ,^'1^ 
we set up the moving TDHFB state vector \(j){q,p, N)) from the static state vector 
\4>{q)) in the following form: 

\4>{q,P,N)) = e^Gfep,7V) = e^PQ('')+-^W|</,(g)), (2-5) 
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where Q{q) and &{q) are one-body operators. Substituting the moving state vector 
\(j}(q,p, N)) of this form into Eqs. (|2-3j) and (|2-4|1 . and comparing the coefficients of 
the first order in p and n, we get 

{m\[QiQ),P{q)]\^iQ))=i, (2-6) 
{^{q)\[0{q),N]\^{q))=i, (2-7) 

where P{q) is the local shift operator of the collective coordinate q, defined by 

Piq) \Hq)) = \<Piq)) ■ (2-8) 

Other commutation relations such as {(l){q) \ [Qiq)i {(piq)) and {(l)iq) \ P{q)] 10(9)) 
are zero. 

The collective Hamiltonian is also expanded up to the second order in p, 
n{q,p,N)^{4>iq,p,N)\H\(P{q,p,N)) 

= V{q) + \B{q)p' + \{q)n, (2-9) 



where 



V{q)=n{q,p,N%=0,N=N, = {m\H\m). i™) 

^(«) = ?? n.. . =-{m\[[H.miQm\m). (2-ii) 

ap'^ p=0,N=No 

= = {<p{q)\ [H,ze{q)] 10(g)) . (2-12) 

The collective potential V{q) expresses static properties of the collective path, while 
the mass function B{q) represents the dynamical property, i.e., inertia of the collec- 
tive motion. The quantity X(q) is regarded as the locally defined chemical potential. 

The basic equations of the ASCC method are obtained by performing an adia- 
batic expansion of the equation of collective path H2-l() with respect to the collective 
momentum p, and requiring the variations to be zero for each order in p. In the 
zero-th order, we obtain the moving-frame HFB equation: 

6{<P{q)\HM{q)\HQ))=0, (2-13) 

where 

HM{q)=H-X{q)N-^Q{q) (2-14) 

is the moving-frame Hamiltonian. 

The first and second orders of the adiabatic expansion of Eq. H2-l() yield the local 
harmonic equations (also called the moving- frame quasiparticle RPA equation).!^ 
They are composed of the following two equations: 

6 (0((?)| [HM{q), Q{q)] - \B{q)P{q) |0(g)) = 0, (2-15) 

5 {m\ [HM{q), ]p{q)] - C{q)Q{q) - ^^[[H^j{q), ^Q(q)],Q{q)] 

-^N\cP{q))=0, (2-16) 
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where C{q) represents local stiffness defined on the collective path, 

The collective variables {q,p), and the collective Hamiltonian TC{q,p,N) are deter- 
mined by solving the ASCC equations ^TT^ . ^TTE^ . and (f^T^ . 

§3. Application of the ASCC to the multi-0(4) model 

In this section, we present an explicit form of the ASCC equations for the 
multi-0(4) model Hamiltonian. The ASCC equations become a very simple form 
for separable interactions. A numerical algorithm to find a collective path is also 
discussed. 

3.1. Multi-0{4) model 

The multi-0(4) model Hamiltonian has been used to test the validity of various 
kinds of theories of nuclear collective motion.^^''^^''^^''^^' In this paper, we extend 
this model Hamiltonian such that it includes the quadrupole-type pairing interaction 
in addition to the monopole-pairing interaction. We write the model Hamiltonian 
in the following form: 

H = ho- ^Go(iU + iit) - ^G2{B^B + BB^) - ]^xb\ (3-1) 

j 

The first term in the right-hand side of Eq. (|3-1() denotes the single-particle Hamil- 
tonian giving the spherical single-particle energy for each j-shell which possesses 
(2i7j)-fold degeneracy (2l?j = 2j + 1). Other terms represent the residual two-body 
interactions; in due order, the monopole-pairing interaction, the quadrupole-type 
pairing interaction, and the quadrupole-type particle-hole interaction. Their inter- 
action strengths are denoted by Go, G2, and X) respectively. The operators ap- 
pearing in this model Hamiltonian are defined in terms of the nucleon creation and 
annihilation operators (cj^, Cjm) by 

it = ^it, Bt = ^d,.B], N = Y^N„ D = Y^d,D„ (3-2) 



where 



m>0 m>0 



m 



with 



1 |m[ < Qj/2, 
-1 \m\ > 



(3-5) 



Here, the operators A and denote the monopole-pair and the number operators, 
while B and D represent the simplified quadrupole-pair and quadrupole particle-hole 
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operators, respectively. These operators contain the factors djajm which simulate 
the basic property of the quadrupole matrix elements {jm\ r^l20 in a schematic 
way. Although they are not real quadrupole operators, wc call them "quadrupole" for 
brevity. Exact solutions (eigen-energies and eigen-functions) of the multi-0(4) model 
are easily obtained by means of the matrix diagonalization method (see Appendix 
A). 

3.2. Quasiparticle representation 

To solve the ASCC equations, it is convenient to work on the quasiparticle ba- 
sis locally defined with respect to the state \(j){q)) on the collective path. For the 
multi-0(4) model, the Bogoliubov transformation to quasiparticle creation and anni- 
hilation operators, aj(g) and aj(g), satisfying the vacuum condition, ai{q) \4>{q)) = 0, 
is written as 

a_i(g) J \Vi{q) Ui{q) J \ c^i J 

Here, the indices zki represent the set of angular momentum quantum numbers 
{j, ±m). 

Using the quasiparticle bilinear operators, 

Aliq) = al{q)alM (3-7) 
Niiq) = al{q)ai{q) + al,{q)a-i{q) , (3-8) 

the nucleon bilinear operators Ajiq) and iVj(g) are rewritten as 

4(9) = u,iq)v,iq) + nf (g)Aj(g) - vf{q)Aiiq) - Uiiq)viiq)Ni{q), (3-9) 
Mq) = ^vHq) + 2ui{q)vi{q){Al{q) + A,{q)) + {u]{q) - vf{q))Ni{q). (3-10) 

The quasiparticle bilinear operators, A|(g), Aj(g), and iVj(g), satisfy the following 
commutation relations 

Ai(g), Aj,(g)] = 6u'{l - Ni{q)), (3-11) 

iV,(g),Aj,(g)] =2Sii,Aliq). (3-12) 

The particle number Nq, the quadrupole deformation D{q), the monopole-pairing 
gap Aoiq), and the quadrupole-pairing gap A2{q) are given by the expectation values 
with respect to the mean-field state vector |</>(g)): 



No = {ct>iq)\N\cf>{q)) = 2^i;2(g), (3.13) 

j>0 

D{q) = {<p{q)\D\4>iq)) = 2Y,dmvUq), (3-14) 

i>0 

Mq) = Go {Hq)\A^ \Hq)) = GoY^u^iqMq), (3-15) 

A2{q) = G2 \Hq)) =G2J2 dmui{q)viiq). (3-16) 



i>0 



Below, we often omit the g-dependence in expressions, like Ai{q) ^ Ai. It 
should be kept in mind that all these quantities are locally defined with respect to 
the quasiparticle vacuum \(l>{q)) and depend on q. 



Effects of Time- Odd Mean Field on Large Amplitude Collective Dynamics 7 



3.3. ASCC equations for separable interactions 

The ASCC equations can be easily solved when the effective interactions in 
the microscopic Hamiltonian are separable. We can always write such a separable 
Hamiltonian in the following form: 



H = ho-\Yl ^sPMpM + 1 Y: ^sF^~^Fs^~\ (3-17) 

where 



2^ " * 2 

s s 



= (F, ± F])/2 = ±Fi^)^ (3-18) 

The superscripts (=b) indicate the Hermitian or anti-Hermitian characters of the bi- 
linear operator F. The multi-0(4) model Hamiltonian under consideration contains 
three kind of residual interactions. The indices, s =1, 2, and 3, to the operators 
F and the interaction strengths Kg represent the monopole-pairing, the quadrupole- 
pairing and the quadrupole particle-hole interactions, respectively: Fs=i = A, Fs=2 = 
B, Fs=3 = D, Ki= 2Go,K2 = 2G2, and K3 = x- 

For the separable Hamiltonian, it is possible to directly derive the ASCC equa- 
tions from the time-dependent variational principle, 

6{cPit)\i^^-h{t)m))=0, (3-19) 
for the self-consistent mean-field Hamiltonian defined by 

h{t) = ho-Y^ im \m + E ^^^^'^ (mp^-^ im) ■ (3-20) 



s 



Here, the exchange terms are neglected as usual for the separable interactions.'^''^ 
The second and third terms in the right-hand side of h{t) represent the time- 
even and time-odd components, respectively. Note that the expectation values 
{(j){t)\Fs ^ |<^(t)) are pure imaginary, indicating that the third term is odd under 
time reversal. Substituting \(j){t)) = e*^'^^'^-* \4'iQ)) aiid expanding Eq. (|3-19|) up to 
the second order in p, we obtain the ASCC equations for the separable Hamiltonian: 

H4>{q)\hM{q)\4>{q)) = 0, (3-21) 
S im [hMiq),Q{q)] - E ^i^i"^ - -Mq)P{q) m) = 0, (3-22) 

s 

S {Hq)\ [hM{q), ]B{q)P{q)] - ^ /{J^ - B{q)C{q)Q{q) - ^ /^(g)^^ 

s s 

+ ^[F]-), (/i(g) - A(g)iV)A]/^i - fN{q)N \^{q)) = 0, 

s 

(3-23) 

where hM{q) denotes the self-consistent mean-field Hamiltonian in the moving frame, 
defined by 

hdq) = kq) - Kq)N - ^Q{q) (3-24) 
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with 

hiq) = ho-Y^ ^sF^''^ {Hq)\F^^^ \m) ■ (3-25) 

s 

{h{q) — X{q)N)A represents the A| and Aj parts of the operator in the parentheses, 
and 

fQ,Ii<l) = -i^s {m\ [Fs^-\QiQ)] , (3-26) 

= ]B{q)P{q)] |<^(g)) , (3-27) 

/ff (9) = -l^s im [[F^^\{Hq) - mN)A],Q{q)] 1^) , (3-28) 
fN{q)=B{q) — . (3-29) 

Note that all matrix elements are real, so that {4>{q) \ Fg ^ \(t^{q)) = 0. In the above 
equations, the quantities /g ] represent the effects of the time-odd components of 
the mean field, which will bring about important effects discussed in the next section. 

The above ASCC equations can be obtained also by inserting the expression 
(|3-17|) into Eqs. (|2-15|1 and (|2-16|) . The exchange terms of the residual interactions 
should be omittedM^ It means that we adopt the Hartree-Bogoliubov (HB) 
approximation in place of the HFB approximation. 

3.4. The moving-frame HB equation 

The moving-frame HB equation (|3-2H) at a given q determines the static TDHB 
state \(j){q)). If we know the operator Q{q), we can solve it using the gradient 

dV 

method. The quantities X(q) and are regarded as Lagrange multipliers which 

oq 

will be determined by the following two constraints. The first constraint is the 
particle number: 

{m\N\^{q))=No. (3-30) 

This constraint specifies the location in the particle number space. The second one 
is 

{4>{q)\QiQ-SQ)\m) = SQ, (3-31) 

which can be readily derived from the canonical variable conditions, 

{m\QiQ)\<t>{q))=0 (3-32) 

and 

{HQ)\^{q)\m) = -l- (3-33) 

Equation (|3-31|) is the condition for the collective coordinate q that the "distance" 
from — 6q)) to \(j){q)) is equal to 6q. In numerical calculations, this quantity, 6q, 
corresponds to a mesh size of the discretized collective path. Provided that Q{q — Sq) 
and Q{q) are known, we can solve the moving-frame HB equation. Because the 
variation in Eq. H3-21() is taken with respect to arbitrary two-quasiparticle creations, 

6\(t>{q))=ala]\(^{q)), (3-34) 
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the two-quasiparticle terms proportional to A| and Aj in the moving-frame Hamil- 
tonian hM{q) should vanish. 

We solve the moving-frame HB equation with the above two constraints by 
means of an algorithm analogous to the gradient methodP Details of this algorithm 
are summarized in Appendix IHl 

3.5. The local harmonic equations 

In order to obtain the collective path, we need to solve the local harmonic 
equations, (|3-22|) and (|3-23j) . and find the operators, Q{q) and P{q), which determine 
the direction of the collective path in the TDHB space. In solving the local harmonic 
equations, we note that the moving-frame Hamiltonian /im(^) is expressed in terms 
of the quasiparticle bilinear operators, aJ, Aj, and Ni, in the following manner: 



fiMiq) = Viq)+Y,Ei{q)Ni, (3-35) 

B,s 

+ J;fW(.)(a1 + A,) + J^F^WiV,, (3-36) 

i>0 i>0 



F^^^ = {<P{q)\F^^^\<Piq))+Fi'-]+Fi 



P^~^ = Y.n:sii)iAl-A,). (3-37) 

Here, 



i>0 





= -Vi), 








= 2diaiUiVi, 


(3-38) 




1 

~ "2' 




1 

= — -djCTj, 


F^A,hi) 


= 0, 


(3-39) 








= -diGiUiVi 


F^^,l{^ 


= diai{uf - 





(3-40) 



Ei{q) = {u^i - vf)iei - ^d^aiD{q) - X{q)) - 2{Ao{q) + diaiA2{q))uiVi. (3-41) 

These quantities are determined by solving the moving- frame HE equation (|3-21j) . 
For later convenience, we define the following quasiparticle bilinear operator: 

= [4+\(%)-A(g)7VU] = 25]4+)(.)(a1- A,) (3-42) 

with 

={2u^{q)vi{q){e^ - xdi<yiD{q) - \{q)) 

-{Ao{q) + dmA2{q)){ul{q) - vKq))] F^^}{i). (3-43) 
The infinitesimal generators can be written as 

Q{q) = Y.Q^{A\+Ai), (3-44) 

P{q) = iY,n{A\-Ai). (3-45) 

i>0 
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We can write down the matrix elements Qi and Pi in terms of j, /p^\ f^l and 
/at by substituting Eqs. (TT^ and (IH^ into Eqs. (IH^ and (TT^ : 

= (2 J- ^2 E ^1^(^)4^ + (2ij;,)L^2 + ^^/^)' (3-46) 

where 

7Vi = 2ni(g)t;,(<7), (3-48) 

/i>''i,. = /g(g) + /g(g). (3-49) 
= ^/Biqjclq). (3-50) 

Substituting Eqs. and (TT^ into Eqs. (|!r77|l and (|!T^ . we 

also have the following relations, 



/Qj=2^^E<.«^- (3-51) 

i>0 

/IS.s = E + • (3-52) 



i>0 

Note that /g 3 = 0. From the canonical variable condition, the orthogonality be- 
tween the collective and number fluctuation modes is required; 

{^{q)\[N,P{q)Mq)) = 2iY,NiPi = 0. (3-53) 

i>0 

Eliminating Qi and Pi from Eqs. and with the use of Eqs. (jTl^ 

and ()3-47|) . we finally reach a dispersion equation, 

S(a;2) • f = 0, (3-54) 

for the quantity f = f(g) = {fq}, fq}, /p^^^i, /]^j^2' /pR,3' /^v}- ^^^^ ^ = i^ij} 
a 6 X 6 matrix whose elements are given by 

Sii = 4GoS«(Fi-\ F^-)) - 1, 5i2 = AGqS^^\F^^},F^^}), (3-55a) 

Si3 = 4GoS(2)(Fi7,\ 5i4 = 4Go5(2)(Fi-\ fW), (3-55b) 

Si5 = 4GoS(2)(F^-), 5i6 = 4Go5(2)(F^-\ iV), (3-55c) 



'S'21 
'S'23 
'S'25 



4G25(2)(F 



A,2 ' V5 
A,2 ' -^A,! ^ 
A,2 ' -^A,3 )•' 



S22 = ^G2S(^\f['>,F\->) 



524 = 4G2S(2)(f(-) 



526 = 4G2S(2)(F\-\Ar) 



A,2 ' -^A,2 /' 
{-) 



(3-56a) 
(3-56b) 
(3-56c) 



531 = 4Go{5«(4+\Fi-)) (3.57^) 
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•532 


= 4Go{5(i: 




^A,2 1 




(3-57b) 


'S'33 


= 4Go{S^^' 








(3-57c) 


'S'34 


= 4Go{S^^' 




^A,2^ 




(3-57d) 


'S35 


= 4Go{S^^' 


y^A,i 


^ A,3 ' 




(3-57e) 


'S36 


= 4Go{S^^' 




N) + 




(3-57f) 



'S41 


= 4G2{5(^' 


(i?2^^ , 




(3-58a) 


'S42 


= 4G2{5(^' 






(3-58b) 


'S43 


= 4G2{5(^' 






(3-58c) 


'S44 


= 4G2{5(i' 




FW) + 5(2)(i?W,FW)}-l, 


(3-58d) 


'S45 


= AG2{S^^' 


y^A,2 ' 




(3-58e) 


'S46 


= AG2{S^^' 


(.^A,2 ' 


iV) + s(2)(i?W,iV)}, 


(3-58f) 



S51 


= 2x{S^'^ 


(4^^ 


F^-l) + u.^Si'\Fi^lFi-i)}, 


(3-59a) 


S52 


= 2x{S^' 






(3-59b) 


S53 


= 2x{S^' 


\^A,3 




(3-59c) 


S54: 


= 2x{S^' 


\^A,3 




(3-59d) 


S55 


= 2x{S^' 


\^A,3 


Fi^)) + 5(^)(i?W,FS5)}-l, 


(3-59e) 


S56 


= 2x{S^' 


v-^A,3 


iV) + 5(2)(4+),Ar)}, 


(3-59f) 



^61 = iu^S^^^ {N, F^^} ) , 562 = a;'5(2) (iV, F^y^^ ) , (3-60a) 

563 = S^^\n, 564 = 5(1) (iV, (3-60b) 

565 = S'^^Hn, 566 = S'^^Hn, N). (3-60c) 

The quantities 5^^^ and 5^^^ are defined by 



i>0 



i>0 



The unknown quantities in the dispersion equation (|3-54p are f(g) and ui'^{q). The 
squared frequency uj'^{q) can be determined by the condition that the matrix S{uj'^{q)) 
cannot have its inverse, 

detS(a;2(g)) = 0. (3-63) 
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When there are many solutions w^(q') satisfying this equation, we choose the smallest 
uj'^{q) solution (including negative value) as the collective mode. Once the io'^{q) 
value and, consequently, the matrix S{q) is specified, the direction of the vector 
f(g) is found. On the other hand, its absolute value is fixed by the normalization 
condition for the collective mode, i.e., 

= 2i^Q,(g)P,((7) =z. (3-64) 

Let us note that we can take an arbitrary scale for the collective coordinate q. This 
means that we can set B{q) = 1 on the collective path without loss of generality. We 
adopt this choice in this paper. Then, uj'^{q) represents nothing but the curvature 
of the collective potential: 

^ d^Vjq) 

There is still an arbitrary choice for the sign of Q{q) and P{q)- This sign specifies 
the "rear" and "front" of one-dimensional collective path. 

3.6. Numerical algorithm for solving the ASCC equations 

The infinitesimal generator Q{q) and P{q), which depends on the quasiparticle 
vacuum \(f){q)), is a solution of the local harmonic equations, while the quasiparti- 
cle vacuum \(p(q)), which depends on Q{q), is a solution of the moving-frame HB 
equation. Thus, the set of ASCC equations require self-consistency and an iterative 
solution. In fact, we need a double iteration for each value of q, because the gradient 
method itself for solving the moving-frame HB equation is an iterative procedure. 
The numerical algorithm utilized in the present work is summarized as follows: 
Step 0: HB state (starting point) 

Solve the static HB equation and choose one of the solution. Then, solve the 
quasiparticle RPA equations and select a collective excitation mode which has 
the lowest frequency w(g'o)- This provides the solution of the ASCC equations at 

q = qo- 

Step 1: Initial setting 

Assume that we have solved the ASCC equations at a position q having the self- 
consistent generator Q[q) and the state 

Set Q^'^\q + 5q) = Q{q) as the initial guess for Q{q + 6q), then start the following 
iteration to find the self-consistent solution at q + 5q. Here, the superscript (n) 
indicates the number of iterations. 

Step 2: Moving-frame HB equation 

Using the operator + 5q) (n > 1), solve the moving-frame HB equation 

at q + 6q, 

6 (^(k^'^^Xq + 6q)\H- A(")iV - + 5q) U(")(g + 5q)^ = 0. (3-66) 

Here, the constraints, Eqs. ()3-30|) and ()3-31() . determine the Lagrange multipliers, 
A^'^) and We use the gradient method described in Appendix El to solve Eq. 
()3-66|) . This determines the moving -frame HB state \(t)'''^\q + 5q)) ■ 

Step 3: Local harmonic equation 

Using \(p^"\q + 6q)) with the Lagrange multipliers X{q + 6q) = X^^\ solve the 
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local harmonic equations, (|3-22|) and (|3-23|) . This determines the infinitesimal 
generator Q^'^\q + 6q). 

Step 4- Self- consistency 

Updating the infinitesimal generator Q^"'\q + 6q), go back to Step 2, and repeat 
Steps 2 and 3 until all quantities at q + 5q converge. 

Step 5: Repetition 

Change q into q + 5q and go back to Step 1. 
Carrying out these iterations, Steps 1-5, we obtain the collective path starting 
from the HB equilibrium point toward one direction (q > qo). The collective path 
toward the opposite direction is obtained by changing the sign of 6q and repeating the 
same procedure, Steps 1-5. In this way, we determine the self-consistent collective 
path passing through the HB equilibrium point. 

We should give an additional remark on Step 3. When we solve the local har- 
monic equations, we set fg I (q) = in the iteration procedure to avoid a numerical 

instability problem.^ We have confirmed that the solutions obtained under this 
assumption satisfy the required self-consistency. Quite recently, we have found that 
this prescription can be justified without loss of generality. This recent progress in 
the ASCC method, including this proof, will be reported in another paper.^ 

The solution of the ASCC equations yields the classical collective Hamiltonian 
(for a constant particle number = A'o): 

n{q,p) = ^p^ + Viq). (3-67) 

We then obtain a quantum collective Hamiltonian by canonical quantization: 7i{q, p) 
1 d 

H{q, --;—). Note that, in this quantization step, there is no ambiguity associated 
i oq 

with the ordering of q and p, because the coordinate scale is chosen such that the 
inverse mass function is unity, i.e., B{q) = 1. 

§4. Numerical calculations and discussions 

4.1. Details of numerical calculation 

We solve the ASCC equations following the algorithm in the previous section 
and determine the collective path embedded in the TDHB phase space. In order to 
investigate the effects of the quadrupole-pairing interaction on the large-amplitude 
collective dynamics, we use the same parameters as those used in Ref. .60 J except 
for the quadrupole-pairing strength G2- We consider a system composed of three 
shells with the spherical single-particle energies, e^^ = 0, e^^ = 1-0, e^^ = 3.5, the 
pair degeneracies, Qj^ = 14, Qj^ = 10, Qj^, = 4, and the single-particle quadrupole 
moment, dj^ = 2, dj^ = dj^ = 1 for each shell. One kind of Fermion is considered in 
this model and the number of particles Nq is set to be 28. This value is the half of the 
total number of shell model states. The quadrupole particle-hole interaction strength 
X is fixed to 0.04, while three values (0.20, 0.16, 0.14) are employed for the monopole 
pairing interaction strength Gq. For each case, effects of the the quadrupole-pairing 
interaction are studied varying its strength G2 (0.00, 0.02, 0.04). These parameter 
range are adopted so that the calculated ratios of the monopole and quadrupole 
pairing gaps simulate the result of realistic analysis of them.^^' The HB calculation 
yields a single spherical minimum for the Go = 0.20 case, and two local minima 
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corresponding to the oblate and prolate equilibrium shapes for the Go = 0.16 and 
0.14 cases. For the latter cases, we present below the results obtained by starting 
from the prolate equilibrium point (D > 0). Of course, we obtain the identical 
collective path even if we start the calculation from the oblate equilibrium point 
{D < 0). Note that the multi-0(4) Hamiltonian possesses the "parity" symmetry 
(invariance under the transformation ajm —o'jm, that is, D —D and B — > —B). 
Therefore, all the quantities are even or odd functions of the quadrupole deformation 
parameter D. 

4.2. Collective path and collective potential 

A nice property of the multi-0(4) model is that we can simulate, by changing 
the ratio Gq/x, the phase transition in a finite quantum system from the single- well 
to the double-well potential. This is analogous to nuclear shape phase transition 
and shape coexistence phenomena. The collective potential V and the pairing gaps 
(Z^O) ^2) are displayed in Figs. ^ and HI In the ASCC method, these quantities are 
calculated as functions of the collective coordinate q, but they are easily converted to 
those as functions of D. See the relation between q and D shown in the left panels of 
Fig. 121 We see that the collective potential change from the single well to the double 
well as the ratio Gq/x decreases, i.e., as the effect of the monopole-pairing interaction 
is weakened. When the oblate-prolate shape coexistence occurs (Go = 0.14,0.16), 
the collective path calculated from one local minimum passes through the other local 
minimum. It is seen that Aq[D) decreases while |Z\2(-D)| increases as D increases. 
Both gaps tend to vanish near the limits, -Dmin = —42 and -Dmax = 42. This is due to 
the small model space and is absent in realistic situations. The quadrupole-pairing 
gap |Z\2(-D)| takes the maximum value at the prolate and oblate HB equilibrium 
points. Because of energy gain associated with the quadrupole pairing, the oblate 
and prolate local minima in the collective potential V{D) become deeper as G2 
increases. 

Significant effects of the quadrupole-pairing interaction are obviously seen in the 
solutions of the local harmonic equations. The squared frequencies represent- 
ing the curvature of the collective potential, are shown in the right column of Fig. |21 
while the sums, Ylii>o\Qi{l)\^ Yli>oWi{l)\^ 1 of ^'^^ quasiparticle compo- 
nents, Qi{q) and Pi{q), are displayed in Fig. |^ They are again plotted as functions 
of D. It is clearly seen that absolute magnitudes of Qi{q) significantly increase, 
while those of Pi{q) decrease with increasing G2. (To avoid complication, their sums 
rather than individual values are presented in Fig. ^) These changes in the micro- 
scopic structure of Q{q) and P{q) lead to the increase of the derivative dq/dD with 
increasing G2 (Fig.|21). Then, this will result in a significant enhancement of the col- 
lective mass M{D{q)) with respect to the coordinate of the quadrupole deformation 
parameter D. These interesting properties of the collective mass are a main subject 
of this paper and will be discussed in the next subsection. 

We have also carried out a constrained HB (CHB) calculations choosing the 
quadrupole operator D as the constraint. Calculated values of V{D), Ao{D) and 
A2{D) were almost indistinguishable from those obtained by the ASCC method dis- 
played in Figs. Hand El Thus, concerning these static mean-field quantities, we see 
no difference between the ASCC and the CHB calculations. It should be noted, how- 
ever, that such a good agreement is due to a simplicity of the multi-0(4) model, i.e., 
it contains only one degree of freedom, D, relevant to the large-amplitude collective 
motion. In reality, many degrees of freedom (a variety of particle-hole excitations as- 
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sociated with shell structures, triaxial deformations, various multipolarities, low and 
high frequency excitations, etc.) would be interwoven to generate large-amplitude 
collective motions. In fact, it has been shown that the self-consistently determined 
collective coordinate operators for low-frequency quadrupole-type collective vibra- 
tions are significantly different from the quadrupole operators (see, for instance, 
Ref. l3n|l). Let us now turn to a discussion on collective mass where we can clearly 
see the merit of the ASCC method. 

4.3. Collective mass 

As mentioned in subsection 13.51 we set the scale of the collective coordinate q 
so as to make the collective mass unity. In order to compare the collective mass 
obtained by the ASCC method with the conventional cranking mass, let us evaluate 
the collective mass as a function of the quadrupole deformation D. This quantity, 
M{D{q)), is readily obtained by transforming the collective kinetic energy as a 
function of the velocity D; 

lB{q)p' = = ^ ^M{D{q))D^ (4-1) 

Mp(g))=('^] ={4Y,diamv^P^{q))-^. (4-2) 

As is well known, the Inglis-Belyaev cranking mass is derived by means of the 
adiabatic perturbation theory^-' and written as 

M.AD)^2^ ^„,^f°^„(„) . (4.3) 

Here, |(/)o(-D)) and Eq{D), respectively, represent the ground state and its energy at 
the deformation D, while \(j)n{D)) and En{D) the (two-quasiparticle) excited states 
and their energies. We can use this formula either treating the deformation D as 
a phenomenological parameter or a self-consistently determined quantity. If the 
ground states are calculated at every points of D by means of the CHB method, 

S {MD)\H - X{D)N - fi{D)D \MD)) = 0, (4-4) 

with the self-consistency conditions for the particle number and the quadrupole 
deformation, 

{MD)\N\MD)) =No, {MD)\D\MD)) = D, (4-5) 
then the following explicit expression for Mct{D) is valid: 



Mer(D) = 2^ 



mm' 

(4-6) 

Hereafter, we call this "CHB-cranking mass." Note that we used a slightly different 
definition of the cranking mass in Ref. 60), in which the CHB self-consistency is 
ignored. 

In Fig. the ASCC mass M{D{q)) are shown as a function of D and com- 
pared with the CHB-cranking mass Mcr{D) for various combinations of the pairing- 
interaction strengths. Go and G2- One may notice that they diverge near Dmin = 
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—42 and Dmax = 42. This behavior indicates that, in the multi-0(4) model un- 
der consideration, it becomes harder and harder to increase D as approaching its 
hmit. Let us then focus our attention on the middle region of D. This is the re- 
gion important for the quantum-mechanical tunneling motion through the barrier 
(small Gq/x case) and for the vibrational motion about the spherical equilibrium 
(large Gq/x case). It is seen that, when the quadrupole-pairing interaction is absent 
(G2 = 0), the magnitudes of the ASCC mass are almost the same as those of the 
CHB-cranking mass. They also exhibit similar deformation dependence. When the 
quadrupole-pairing interaction is switched on, however, a significant difference be- 
tween the ASCC and the CHB-cranking mass appears: The ASCC mass significantly 
increases as G2 increases, while the CHB-cranking mass hardly changes. 

The origin of this different behavior between M{D{q)) and Mcr(-D) may be 
understood in the following way: The CHB-cranking mass is derived by neglecting 
the contribution of the residual interaction to the collective mass, so that effects of 
the quadrupole-pairing interaction are taken into account only in the static quantities 
like the quadrupole-pairing gap A2{D). On the other hand, the effects of the residual 
interaction on dynamics are taken into account in the ASCC method through the 
time-odd components of the mean-field that change sign under the time reversal p 
—p. The time-odd part of the mean-field Hamiltonian h{t) represents the change in 
the self-consistent mean-field associated with the dynamical motion. It is well known 
that this part is indispensable for reproducing the correct center of mass for the 
translational motions.^' On the other hand, it is also known that both the monopole- 
pairing and the quadrupole particle-hole interactions do not contribute to the time- 
odd part.l^'-^^''^^' Thus, only the quadrupole-pairing interaction contributes to it in 
the case of the multi-0(4) model Hamiltonian under consideration. 

The remarkable difference between the ASCC mass and the CHB-cranking mass 
displayed in Fig. [S] clearly indicates the importance of the quadrupole-pairing in- 
teraction on collective dynamics. This shows a striking contrast to its effect on 
static properties. Namely, the major properties of the collective potential energy 
curve V{D) are determined by the competition between the monopole-pairing and 
quadrupole particle-hole correlations. There, A2{D) is much smaller than Aq{D), 
and the quadrupole pairing plays only a minor role. However, it turns out to play a 
major role for dynamical properties of collective motion. 

4.4. Excitation spectra and transition matrix elements 

We calculate the eigen-energies and wave functions ^kio) for the quantized 
large-amplitude collective motion by solving the collective Schrodinger equation: 



and the boundary conditions, "^kilmm) = "^kilmax) = 0, where q^m and gmax are the 
minimum and maximum values of q along the collective path obtained (see Fig. 3). 
The quadrupole transition matrix elements are evaluated by 




(4-7) 



with the orthonormalization 




(4-8) 




(4.9) 
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where the deformation D{q) on the collective path is defined by Eq. (|3-14|) . 

Figure IHl displays the results of the ASCC calculation for excitation spectra 
and quadrupole transition matrix elements between low-lying states. For Gq = 
0.20 we obtain anharmonic vibrational spectra about the spherical equilibrium. In 
contrast, the collective potential V{D) for Go = 0.14 and 0.16 possesses two local 
minima corresponding to the oblate and prolate shapes, and the spherical point 
(D = 0) becomes the barrier top. In the Go = 0.14 case, this barrier is high. 
Consequently, a ground state doublet similar to the well-known parity doublet in 
the double well potential appears. In the present multi-0(4) model, the doublet 
corresponds to the symmetric and anti-symmetric superpositions of the oblate and 
prolate ground states, and its energy splitting provides a sensitive measure of the 
quantum tunneling effect through the potential barrier. On the other hand, in the 
Go = 0.16 case, the barrier is rather low so that the spectrum exhibits a transient 
feature toward the doublet pattern mentioned above. In the quantum spectra of 
Fig. El we can clearly identify the effects of the increase of the collective mass due 
to the quadrupole pairing. First, the vibrational excitation energy decreases as G2 
increases. Second, the energy splitting of the doublet decreases with increasing G2 
indicating that the tunneling motion becomes harder as the collective mass increases. 
Thus, for the combination of Go = 0.14 and G2 = 0.04, we obtain a doublet of 
excited states in addition to the ground-state doublet. The wave functions of the 
excited-state doublet as well as the ground-state doublet are displayed in Fig. 9. 
This figure clearly indicates that the excited doublet corresponds to the symmetric 
and anti-symmetric linear combinations of vibrational excitations about the oblate 
and prolate local minima. We can also confirm that the amplitude in the barrier 
region indeed decreases with increasing G2. 

In Fig. [3 the results of the exact matrix diagonalization of the microscopic 
multi-0(4) Hamiltonian are presented. For every combinations of the interaction 
strengths, Gq,G2 and X) the excitation spectra and the transition matrix elements 
obtained by the ASCC method agree in a very good approximation with the results 
of the exact calculation. It is remarkable that the ASCC calculation succeeds in 
describing the gradual change of the quantum spectra associated with the phase 
transition of the finite system from the spherical shape to the oblate-prolate shape 
coexistence. 

Let us make a severe quantitative comparison concerning the energy splitting of 
the ground-state doublet in the Go = 0.14 case. The splittings obtained in the ASCC 
method are 0.043, 0.012, and 5x10"^ for G2 = 0.00,0.02, and 0.04, respectively. 
The corresponding values obtained by the exact diagonalization are 0.091, 0.020, 
and 3xl0~^. It should be noted here that the energy splitting under discussion 
is a very small quantity associated with the barrier penetration for which even a 
slight error in the collective mass will result in an error of the order of magnitude. 
Therefore, the agreement within the factor of two indicates that the collective mass 
evaluated by the ASCC method is very reliable. It should also be emphasized that 
this tunneling motion is a large-amplitude collective motion associated with the 
major rearrangement of microscopic configurations of many particles and that the 
collective mass represents the inertia of this motion of many-body system as a whole. 
Thus, the correct evaluation of the collective mass is a highly non-trivial task. 

Let us now investigate how the difference between the ASCC mass and the CHB- 
cranking mass discussed in the previous subsection affects the excitation spectra and 
transition matrix elements. Adopting the Pauli quantization prescription, we obtain 
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the collective Schrodinger equation 




(4-10) 



for the wave functions (D) that incorporate the metric factor Mct{D)^^^ such 
that the orthonormalizations are given bjJ^ 




^r*{D)^riD)dD = 6ki. (4-11) 



We solve this Schrodinger equation under the boundary conditions, 5?'^'^'^ (-Dmin) = 
!Z'^^'^^(-Dmax) = 0. The quadrupole transition matrix elements are evaluated by 



As mentioned in subsection 14. 21 the ASCC collective potential V{q) almost coincides 
with the potential Vchb(-C) calculated by means of the CHB method, so that the 
major difference in the quantum spectra can be attributed to the difference in the 
collective mass. The results are displayed in Fig. |HJ When the quadrupole pairing 
is absent {G2 = 0), we see rather good agreements with the ASCC and with the 
exact solutions. However, the excitation spectra and transition matrix elements are 
almost unchanged when the quadrupole-pairing interaction is switched on and G2 
increases. This is very different from the ASCC results and the exact solutions. We 
can confirm this point also in the wave functions displayed in Fig. |H1 The amplitudes 
in the barrier region change only little even when G2 increases, in contrast to the 
ASCC wave functions. It is obvious that this failure in taking into account the effect 
of the quadrupole pairing originates from the fact that the CHB-cranking procedure 
ignores the time-odd mean- field contribution to the collective mass. 



The multi-0(4) model is a simple model to simulate the phase transitions in 
finite quantum systems from the spherical shape to the oblate-prolate shape coexis- 
tence. We have applied the ASCC method to this model and studied the collective 
mass of the many-body tunneling motion through the potential barrier between the 
oblate and prolate local minima. Comparing with the results of exact diagonal- 
ization, we have shown that the ASCC method succeeds in describing the gradual 
change of the excitation spectra from the anharmonic vibration about the spherical 
equilibrium to the doublet pattern associated with the deformed double-well poten- 
tial possessing the oblate-prolate symmetry. It was found that the collective mass 
significantly increases due to the quadrupole-pairing contribution to the time-odd 
component of the moving mean field. We have also shown that the CHB-cranking 
procedure underestimates the collective mass, because the contribution from the 
time-odd component is disregarded there. 

Along the approach developed in this paper, we shall evaluate in a forthcoming 
papeJ^ the contribution of the time-odd components to the collective mass of large 
amplitude collective motion associated with the oblate-prolate shape coexistence in 
the ^^Se and ^^Kr region. 




(4-12) 



§5. Conclusions 
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Appendix A 

Exact diagonalization of the multi-0{4) model Hamiltonian 

It is possible to construct the following two sets of SU{2) generators for each 
j-shell from the basic operators of the multi-C>(4) model: 













(A-1) 








= 




(A-2) 












(A-3) 
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The operators {Kj^, Kj-, Kjq} and {Lj^, Lj^, Ljq} satisfy the commutation rela- 
tions of the SU (2) algebra: 

[K+,K-] = 2K^6,y, [L+, LJ,] = 2L%y, (A-4) 

[K^, Kp] = ±Kf6jj,, [Ll Lf,] = ±Lf6j,,. (A-S) 

All commutation relations between K^s and L's are zero. Thus, we can construct 
the orthogonal basis vectors of the model space as 

\nK,nL) ^ll\^K„nL^) ^ll{K+r-.{L+r.\0), (A-G) 

j j 

where {nXj^nLj) are integer numbers satisfying < nxpUi^ < f2j/2 and Ylji^Kj + 
nL,) = NQ/2. 

In terms of these SU{2) generators, the multi-0(4) Hamiltonian is expressed as 
H = Y^ eO{2(i^o + L") + n,}-2xY^ d,d,{Kf - L^){K^ - L") 

i ^(Go + G2d,d,)iK+Kr + K+Kr + L+LJ + LrLp 

i ^(Go - G2d,d,){KrLj + K+LJ + LrK+ + L+Kr). (A-?) 



2 



The operators, Kp,K- and K^, act on the basis |ni^^. ,nL^) as 



|?^ii-j,raL,) = J {uK, + 1) f - J^x, ) l^-i^, + l,nL,) , (A-8b) 



Similar equations hold for L'j',LJ, and L^. The matrix elements of the multi-0(4) 
Hamiltonian, n£,'| if [nx, «,£,), can be calculated by using Eqs. (|A-7|) and 

HA-8|) . Diagonalizing this matrix, we obtain exact eigen-energies and eigen-states. 
For the parameters given in the text, the dimension of this Hamiltonian matrix is 
1894. The quadrupole transition matrix elements between eigen-states, \(j)a) and 
\(pi3), are given by 

{cP^\D\cPfs)= Yl '^dA^K,-nL^)C^*j^^riAK,nL^ (A-9) 
riKriL 

where .j^^ are expansion coefficients in the SU{2) basis, 

\<t^a)= C^K,nL\^K,nL) . (A-10) 

riKfiL 
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Appendix B 

Gradient method for the moving-frame HB equation 



We solve the variational equation of the following form using the gradient 
method:''^ 



5{(t>{q)\H-\N -^JiQ\m)=^. 



with the constraint conditions for the number operator N and a one-body operators 
R, 



{m\ N m)) = No {<l>{q)\ R \<j)iq)) = Ro- 



(B-2) 



Here, R is an arbitrary one-body operator which, in general, may be different from 
Q. 

Let \4>^^\q)) be the state vector at the iterative step k. Using the quasiparti- 
cle bilinear operators, A|(g) and Ai{q) which satisfy Ai{q) \ (t>^'^\q)) = 0, we then 
generate the state vector at the {k + l)-th step in a form of unitary transform of 
{q)) as follows: 

(B-3) 



(B-4) 



with the anti-Hermitian operator 

Z{q) = Y,Zf\q){Al{q)-Mq)). 



i>0 



It should be noted that the normalization is preserved during the iteration. As- 
suming that Z(q) is small, we expand the energy difference between \4>^''^^\q)) and 
\(j)^''\q)) as follows: 



AE 



H-XN -fxQ 



H-XN -ijQ 



<P^>^Hq) [H-XN-iJiQ,Z{q)] (p^^\q)) + 0{Z^) 
J2{Hf\q) - AAf (,) - ^^Qr{q))Z!%) + 0{Z'). 



i>0 



If Z^ (q) is chosen as 



Zf'iq) = -ATiHfiq) - XN^iq) - f^Qf{q)) 



r20 



20/ 



-)20/ 



(B-5) 
(B-6) 
(B-7) 

(B-8) 



with a positive step size AT, AE is negative in each iteration. The constraint 
conditions ()B-2|) can also be expanded up to the first order in Z(q), 



-1) 



R 



kik+l) 



cl^C^Hq) R + [R,Z{q)] c/^^'Hq)^ 

+ 0{Z^), (B-9a) 
N \<P^'+'\q)) = {^^"Hq)] N l^^'Kq)) + {^^"Kq)] [A, Z{q)] \<P^''Kq)) 

+ 0{Z^). (B-9b) 
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Substituting Eq. ()B-8|) into Eq. ()B-9|) . we obtain the following equation which de- 
termine the Lagrange multipliers A and /u: 

i i 
_ i i 

where N^^, Q?^ and Hf^ denote the coefficients of the two-quasiparticle creation 
and annihilation parts of the operators N, R, Q and H, respectively, while the quan- 
tities R and N represent {4>^^\q) \ R \(t)^'^\q)) and {(t)^^\q) \ N \(l3^^\q)), respectively. 

When solving the moving-frame HFB equation H3-21() with the constraint con- 
ditions, (|3-3flj) and (|3-31|1 . we use Q{q — 5q) for R{q). This operator is defined at 
q — 5q m terms of the quasiparticle bilinear operators, A\{q — 6q) and Ai{q — Sq), 
which satisfy Ai{q — 6q) — 5q)) = 0. Thus, at each iteration step, it is necessary 
to rewrite Q{q — 6q) in terms of A|(g), Ai{q) and Ni{q), defined with respect to 
I (g)> at q. 

R{q) = Qiq - 5?) = E ^^''^^ - - +Mq- Sq)) (B-11) 

= R'^q) + RlHQ)N^{q) + i?f (g)(Aj(<?) + A{q)). {B-12) 

i>0 

The explicit expression for the coefficients Rf^ (g) , which we need in solving Eq. ()B-10|) , 
is 

Rl\q) = Qfiq - 5q) [{uM - 5q)^ -V,{q- 5q)^){uj{q) - vfiq)) 

+4:Ui{q - 5q)vi{q - 5q)ui{q)vi{q)] . (B-13) 



A 



Rq — R 
2AT 

No-N 
2AT 



20^20 



(B-10) 



The above procedure is repeated until convergence is attained. In fact, this imple- 
ments Step 2 in the iterative algorithm described in subsection 13.61 
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g. 2. Monopole-pairing gaps Aq {left column) and quadrupole-pairing gaps A2 {right col- 
umn), plotted as functions of D. The upper, middle and lower rows display the results for 
Go = 0.14,0.16, and 0.20, respectively. In each panel, the results for G2 = 0.00,0.02,0.04 are 
compared. 
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Fig. 3. Left column: Relation between tlie collective coordinate q and the quadrupole deformation 
D{q) = {'Pi<l)\ D \4>{q)) . The point q = corresponds to the HB equilibrium, which is the 
starting point of the numerical calculation. Right column: Squared frequencies co'^{q) of the 
local harmonic equation, plotted as a function of D. Note that they arc negative, i.e., uj{q) is 
pure imaginary, in the region where the curvature of the collective potential is negative. The 
upper, middle and lower rows display the results for Go = 0.14, 0.16, and 0.20, respectively. In 
each panel, the results for G2 = 0.00, 0.02, 0.04 are compared. 
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Fig. 4. The sums, X]i>o \Qiil)\^ ^^'^ X]i>o of the two quasi-particle components, Qi{q) 

and Pi{q), of the infinitesimal generators, Q{q) and P{q), plotted as functions of the quadrupole 
deformation D. The upper, middle and lower rows display the results for Go = 0.14, 0.16, and 
0.20, respectively. In each panel, the results for G2 = 0.00, 0.02, 0.04 are compared. 
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Fig. 5. The ASCC collective mass M{D{q)) (left column) and the CHB-cranking mass Mcr{D) 
{right column) as functions of deformation D. The upper, middle and lower rows display 
the results for Go = 0.14,0.16, and 0.20, respectively. In each panel, the results for G2 = 
0.00, 0.02, 0.04 are compared. 
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Fig. 6. Excitation spectra calculated with the ASCC method. The upper, middle and lower rows 
display the results for Go = 0.14, 0.16, and 0.20, respectively. In each row, the results for G2 = 
0.00, 0.02, 0.04 are compared. The numbers adjacent to vertical lines indicate absolute values of 
the transition matrix elements. The matrix elements between the doublets are indicated beside 
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Fig. 7. Excitation spectra calculated with the exact diagonalization. See the caption of Fig. |S1 
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Fig. 9. Wave functions of eigenstates of the collective Schrodinger equation, Eq. H4-l()|l . for Go = 
0.14. Those obtained with the ASCC (ClfB-cranking) method are shown in the left (right). In 
order to make a comparison, they are plotted as functions of a common parameter D. Namely, 
those of the ASCC method are defined by 9k{D{q)) = ^k{q)^/dq/dD = Wk{q)M{D{q)Y/^ and 
normalized as J \'^k{D)'^dD — 1. The first, second, third and fourth rows display those of the 
ground state, the first, second, and third excited states, respectively. In each panel, the results 
for G2 = 0.00, 0.02, 0.04 are compared. 



